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Abstract 

We extend the semiclassical study of the Neumann model down to the 
deep quantum regime. A detailed study of connection formulae at the 
turning points allows to get good matching with the exact results for the 
whole range of parameters. 
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1 Introduction 



The Neumann model is an interesting example of integrable model, both clas- 
sically and quantum mechanically, which describes the motion of a point on a 
sphere subject to a harmonic potential. In the quantum lot of effort 

has been devoted to study the semi-classical approximation [T|, since it deals 
with algebro-geometric objects which appear naturally in the classical theory. 
However this study becomes more interesting and concrete when one has pre- 
cise numbers to state effectively the Bohr-Sommerfeld conditions and compare 
them with exact spectra. In a previous works [3 E], we showed how to solve 
numerically the separated Schrodinger equation to obtain numerical values of the 
energy and compared this to a semiclassical computation valid for a large radius 
of the sphere. This case corresponds to the localization of the particle around the 
antipodal minima of the potential. These localized states come in pairs which 
are split by the tunneling probability between the two poles. The two states are 
either symmetric or antisymmetric under exchange of the poles. 

When the sphere radius shrinks, the states become delocalized and the de- 
generacy is lifted by a finite quantity. In the limit of zero radius, the potential 
becomes irrelevant and the energy is the one of the free particle on the sphere 

1 / 2j ( j + 1). The symmetric and antisymmetric solutions go to states with j dif- 
fering by one. When flowing between the large radius and small radius cases, the 
matching conditions used in the WKB analysis must change, in particular to ex- 
plain this degeneracy lifting. We study in this work the evolution of the matching 
phase which either goes from ir/A to through ir/8 for a symmetric wave function 
of from 7r/4 to it/ 2 through 3n/8 for an antisymmetric wave function. 

In a first part, we study the zero radius limiting case, where the energy is 
reproduced by the WKB method up to a small constant. Then we detail the 
numerical evolution of the phase. In a last part, we explain this evolution through 
the asymptotic analysis of a simplified model. 

2 Semiclassical analysis of spheroidal harmon- 
ics. 

Even if the Neumann model has been generalized to arbitrary dimensions, there 
are no really new phenomena appearing above the three dimensions to which 
we limit ourselves in this work. The parameters of the Neumann model are the 
oscillator strengths which can be reduced without loss of generality to be 0, 1 and 
y. The Schrodinger equation separates using the Neumann coordinates, yielding 
the one-dimensional equation: 
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In this equation, f 2 and ^3 are the eigenvalues of the conserved quantities which 
have to be determined and v is r 2 / (Ah 2 ) where r is the radius of the sphere. The 
energy is E = 2ft 2 (/ 2 + yfz)- In the sequel, we take ft — 1. The points 0, 1 and 
y are regular singularities of eq. with exponents and 1/2, corresponding to 
solutions with monodromies ±1. In order to recover a well defined solution on 
the sphere with definite parity properties under the three possible reflexions, the 
solution must have definite monodromies at the three singularities [2|. An even 
solution is thus a solution with monodromy 1 at 0, 1 and y, that is a function 
which is analytic on the whole complex plane. If we want an odd solution under 
x\ — > —Xi, we factor \/t in \1/ and search for an analytic solution to the equation: 



/ d 2 1/3 1 1 \d 

\~ctt 2 + 2 \t + t - 1 + t-y) dt ~ 

-f2-fs + 1/4 + !/(%) h ~ 1/4 / 3 -V(4y) \ 
t t-1 t-y J 



*(t) = (2) 



Similar equations can be written for the six other possible combinations of pari- 
ties. 

In our previous work [3], we showed how the semiclassical analysis nicely fits 
the exact numerical spectra of the Neumann model in the large v limit, when the 
point is confined around the poles. To extend the analysis to the whole range 
of values for v, we first consider the v = case, corresponding to the spheroidal 
harmonics. The energy is known exactly in this case, since it is simply the energy 
of a free particle on the sphere, \j(j + 1). 

In this case, the potential simplifies since the numerator is of degree one: 

v _ _E t-b 

~~2t(t-l)(t-y) U 

The Bohr-Sommerfeld integrals are thus elliptic integrals, with the four special 
points 0, 1, b and y corresponding to the vertices of the period rectangle, and a 
pole at infinity. We have two cases to study according to the position of b with 
respect to 1. 

The quantification conditions, in the case with all monodromies equal to 1, 
are either ones of: 

b>l, [ Wdt = mir, [ V Wdt = {n + -)7i (4) 
Jo Jb 4 

6<1, / Wdt = (m+^)vr, f V Wdt = nn (5) 
Jo 4 J l 

These quantification conditions involve the usual 7r/4 phase factor coming from 
the asymptotics of the Airy function around the turning point b, see e.g. j3]. 

Adding the two Bohr-Sommerfeld integrals, one gets a condition which does 
not depend on the position of b with respect to 1. In fact, this sum can be 
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computed and does not depend on b. Indeed, the integrand is an abelian integral 
of the third type with simple poles at the two points over t — oo. The real 
part of the integral on the real axis is just the sum we want to evaluate, since 
the integrand is purely imaginary on the complementary segments. Hence the 
residue theorem allows to calculate this sum. The residue is i\J E/2 and since 
the integration contour runs through the pole, we get a contribution of m times 
the residue. The final quantification condition is therefore: 

n^/Wj2 = (m + n + ^)tt (6) 

The semiclassical energy thus depends only on the total number of excitations and 
not on the individual values of m and n. The angular momentum is j = 2(m + n) 
so that E = l/2(j + l/2) 2 = l/2j(j + l) + l/8. This differs by merely 1/8 from the 
exact result and is exactly the value of L 2 which must be used in a semiclassical 
treatment of spherically symmetric potential according to Langer 

When b varies, we should be able to pass smoothly from the case of eq. (JH) 
to the one of eq. (jSj). The additional phases therefore go from to 7i~/4 and from 
7r/4 to 0, but their sum should remain constant for the semiclassical energy to be 
completely independent of b. In particular, in the case b = 1, we should have for 
symmetry reason the same it/ 8 phase on both sides. We shall show below that 
this is indeed the case. 



3 Numerical investigations. 

In order to study the variation of the phases appearing in the connection formulae 
at the turning points, we numerically solve the equation using the methods 
of [2]. From the values of f 2 and fs, we calculate the zeros b\ and 62 of the 
potential, which can be written: 

= _ («->.)(«->.) (7) 

t(t - l)(t - y) ' 

The Bohr-Sommerfeld integrals are taken on intervals where the potential is 
positive, with boundary points in the set {0, 1, y, bi, b 2 }. We numerically evaluate 
them and plot their differences with naive expectations in units of n. In the 
following plots, we study the case where the excitation numbers are n = in the 
interval [0, 1] and m = 4 in the interval [l,y = 2]. We obtain similar results for 
other values of the parameters n, m and y. The plot (Figure P) clearly shows that 
the additional phase in the interval [1, 2] goes from when 6 2 is smaller than 1 
to a value approaching 7r/4 in the other case, with the variation stronger when 62 
crosses 1. To isolate the dependence on bi, we plot the sum of the two integrals 
which should not depend on the position of 62 with respect to 1. It varies from 
7r/4 for negative b\ to n/2 for positive b\. 
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Figure 3: Variation of bi, &2 and phases as function of v. 

In order to show the details of the transition region, we next plot (Figure EI) 
the relevant phases as functions of b\ or b 2 — 1. We see that the dependence of 
the phase is not monotonic, but there is a small "glitch" just at b\ — or b% — 1. 
We shall see below that this can indeed be explained by analytical studies. We 
also see that the value of the phase is not precisely it/ 8 at this point of transition. 

We further study with different monodromies and values n — 1 and 

m = 3 for the number of zeros in each of the subintervals. Precisely in our 
example we choose monodromy -1 at and 1 and monodromy +1 at y=2, and 
obtain the plot (Figure IHJ). In this case, the additional phases decrease from n/2 
to 7r/4 through 37r/8 in the transition region. This difference of behavior with the 
case of monodromies +1 perfectly explains the lifting of degeneracy which occurs 
at small v between the two degenerate levels at large v. The relative plot (Figure 
@| shows that in this case the transition is smoother, which will be explained in 
the next section. 

4 Theoretical study of the transition region. 

To understand the behavior of the phases in the transition region we consider for 
example the region around the singularity t=0 and formulate a simplified version 
of eq.fpP) for the case of monodromy 1 or of eq.(|2*|) for the case of monodromy 
-1. The simplified model uses that the terms having poles at t=l or t=y will 
not vary much as long as t remains close to 0, and doesn't approach the next 
pole 1. Hence we simply replace these polar terms by constants, and study the 
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equations, where fi = v — f 2 — f^. 

£ + (s^)5-(4 + *))'W" < 8 > 

Equations of this form can be readily related to the confluent hypergeometric 
differential equation, since there is clearly a regular singularity at t — and an 
irregular one at t — oo. 

Recall the differential equation satisfied by the confluent hypergeometric func- 
tions 

d 2 w . efau , . 

+ = ° (10) 

The solution analytic at the regular singularity z = is: 

M( ,M) = i + T + ^-L + ... 

One can find its asymptotics at the irregular singularity z = oo by relating it 
to the Whittaker functions which have simple asymptotics there, through the 
Mellin-Barnes transformation j^]. The result is given [7j by; 

M(a,b,z) e ±i™ z ~a T.\ e 2 ^" 6 



r(6) " r(6 



fs( i+o y) + w( i+ °y) (n) 



where the plus sign is taken in the sector —tt/2 < arg^ < 3ir/2 and the minus 
sign in the sector — 3n/2 < axgz < — ir/2. 

To get eqs. (|8|9jl into the form eq. ffT0|) we take ^/(t) = exj>(at)w(z) and t = (3z 
and observe one gets exactly the confluent hypergeometric equation taking: 



A A 2 
a = -— + i\B-—, f3 



2 V 4 - A 2 

Then one obtains, starting from eq. (jHJ) that: 

a = \-P(fi-\A), b= 1 - 
and similarly starting from eq. (JOJ) that: 

Finally remembering that we need precisely the solution which is regular at t=0 
in both cases, the corresponding solution is thus M(a,b, z), whose asymptotic 
expansion we know from eq.Qlip. and which can be directly compared with the 
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semi-classical approximation. Of course the validity of this procedure rests on the 
hypothesis that t is large enough so that one can apply the asymptotic analysis, 
still small enough that the next pole, t — 1 remains distant. This is the case 
when (3 is small, and is thus justified in our case for large B, that is when fi and 
/3 are large enough. Note that in this case A is small compared to B. In the two 



cases, a takes the form \ — iF or | — %F' . 



Applying the asymptotic formula (fTTj) we get for large positive t, respectively: 

1 7T 1 

tf(t) ~ -^cos((jt- - + Fln(2ut)-%lnr(- + iF)) (12) 

tt(t) ~ -rj-cos{ujt-^- + F'ln{2ujt)-Q\nr(- + tF')) (13) 
t I 8 4 

with co> = \J B — A 2 /A. We then have to compare these results with the phase 
we obtain from the WKB analysis of eqs. (|S]L which covers the cases of the two 
monodromies. The semiclassical action is given by: 



f \ ^^dt' = y/t(Bt + A) + -4= ln(l + 2Bt/h + 2 V(l + Bt/ h){Bt/ f x )) 
Jo ' 1 2v f) 

In the limit of large Bt/ fx, this is simply 

ut + F \n(2ut) + F -F ln(F) 

with the same quantities and F defined above, but taken with A = 0. If A = 0, 
we deduce that the phase difference between the semiclassical result and the exact 
one are therefore: 

A0+ = - + 31nr(- + iF)-Fln(F) + F (14) 
8 4 

3-7T 3 

A0_ = — + 31nr(- + zF') -F'ln(F') + F' (15) 
8 4 

In fact, F = F' when A = 0. The variation of the phases A0± are plotted 
in (Figure EJ), which shows that this analytical study recovers the limit of the 
phases found numerically in the preceding section. This plot shows a striking 
similarity with the above numerical results, complete with the 'glitch' around 
in A0 + , which comes from the singularity of the derivative of Fln(F) at F = 0. 
One could worry on the possibility to make a proper match when A is non zero. 
However, the difference on u is quadratic in A and is further divided by VB, so 
that the corresponding phase difference cannot become significant over the range 
necessary to obtain a proper match. The term proportional to the logarithm of 
out is linear in A, but the slow variation of the logarithm makes it also benign, so 
that the only contribution which is significant is in the phase of the T function. 
This should explain the fact that when the bi cross their critical line, which 



9 



corresponds to fx = 0, we obtain phases which are different from the tt/8 or 37r/8 
which we obtain in the simpler treatment with A = 0. Changing the sign of t 
corresponds to keeping B while changing the signs of f\ and A. Our formulae 
for the A<p± show that the phases on both side of a singularity add up to 7r/4 or 
37r/4 as we supposed in section 2 for the invariance of the semiclassical energy at 
v = 

5 Conclusion. 

The combined results of [3] and the present work show that it is possible to re- 
produce the entire spectrum of the quantum Neumann model with semiclassical 
methods in the whole range of parameters. The accuracy remains good even for 
low lying levels, but finite differences remain between the semiclassical and exact 
results. The new feature in the realm of semiclassical studies is that the known 
facts about the asymptotic behaviour of confluent hypergeometric functions al- 
lows to properly model the transition region between different types of boundary 
conditions. This introduces contributions which cannot be reduced to a discrete 
Maslov index. 
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